Why is Tin so soft? 
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Abstract 

The distribution of isoscalar monopole strength in the neutron- even 112 ~ 124 Sn -isotopes has been 
computed using a relativistic random-phase-approximation approach. The accurately-calibrated 
model used here ("FSUGold") has been successful in reproducing both ground-state observables 
as well as collective excitations — including the giant monopole resonance (GMR) in 90 Zr, 144 Sm, 
and 208 Pb. Yet this same model significantly overestimates the GMR energies in the Sn isotopes. 
It is argued that the question of "Why is Tin so soft?" becomes an important challenge to the 
field and one that should be answered without sacrificing the success already achieved by several 
theoretical models. 
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The compression modulus of nuclear mat- 
ter (also known as the nuclear incompressibil- 
ity) is a fundamental parameter of the equa- 
tion of state that controls small density fluc- 
tuations around the saturation point. While 
existing ground-state observables have accu- 
rately constrained the binding energy per nu- 
cleon (B/A~— 16 MeV) and the baryon den- 
sity (p ~ 0.15 fm -3 ) of symmetric nuclear 
matter at saturation, the extraction of the 
compression modulus (K) requires to probe 
the response of the nuclear system to small 
density fluctuations. It is generally agreed 
that the nuclear compressional modes — par- 
ticularly the isoscalar giant monopole reso- 
nance (GMR) — provide the optimal route to 
the determination of the nuclear incompress- 
ibility [2]. Moreover, the field has attained a 
level of maturity and sophistication that de- 
mands strict standards in doing so. It is now 
demanded that the same microscopic model 
that predicts a particular value for the com- 
pression modulus of infinite nuclear matter 
(an experimentally inaccessible quantity) be 
able to accurately reproduce the experimen- 
tal distribution of monopole strength. 

Earlier attempts at extracting the com- 
pression modulus of symmetric nuclear mat- 
ter relied primarily on the distribution of 
isoscalar monopole strength in 208 Pb - - a 



heavy nucleus with a well developed giant res- 
onance peak 0, H|. However, as was pointed 
out recently in Refs. [3, ll| - - and confirmed 
since then by several other groups 0, Q, [1] 
- the GMR in 208 Pb does not provide a 
clean determination of the compression mod- 
ulus of symmetric nuclear matter. Rather, 
it constraints the nuclear incompressibility of 
neutron-rich matter at the particular value of 
the neutron excess found in 208 Pb, namely, 
b = {N-Z)/A = 0.21. As such, the GMR 
in 208 Pb is sensitive to the density depen- 
dence of the symmetry energy. The sym- 
metry energy represents a penalty levied on 
the system as it departs from the symmetric 
limit of equal number of neutrons and pro- 
tons. As the infinite nuclear system becomes 
neutron rich, the saturation density moves to 
lower densities, the binding energy weakens, 
and the nuclear incompressibility softens [||. 
Thus, the compression modulus of a neutron 
rich system having the same neutron excess 
as 208 Pb is lower than the compression modu- 
lus of symmetric nuclear matter. We note in 
passing that the symmetry energy is to an ex- 
cellent approximation equal to the difference 
between the energy of pure neutron matter 
(with 6=1) and that of symmetric nuclear 
matter (with 6 = 0). 

The alluded sensitivity of the distribu- 
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tion of isoscalar monopole strength to the 
density dependence of the symmetry energy 
proved instrumental in resolving a puzzle in- 
volving K: how can accurately calibrated 
models that reproduce ground state data as 
well as the distribution of monopole strength 
in 208 Pb, predict values for K that differ by 
as much as 25%? This discrepancy is now 
attributed to the poorly determined density 
dependence of the symmetry energy Q]. In- 
deed, models that predict a stiffer symmetry 
energy (one that increases faster with den- 
sity) consistently predict higher compression 
moduli than those with a softer symmetry 
energy. Thus, the success of some models 
in reproducing the GMR in 208 Pb was ac- 
cidental, as it resulted from a combination 
of both a stiff equation of state for symmet- 
ric nuclear matter and a stiff symmetry en- 
ergy [5]. Since then, the large differences 
in the predicted value of K have been rec- 
onciled and a "consensus" has been reached 
that places the value of the incompressibil- 
ity coefficient of symmetric nuclear matter at 
fT = 230±10 MeV @, @, M, M ■ 



An example of how this consensus was 
reached is depicted in Fig. [1] where the dis- 
tribution of isoscalar monopole strength in 
90 Zr, 116 Sn, 144 Sm, and 208 Pb at the small 
momentum transfer of q — 45.5 MeV (or 
g = 0.23 fin -1 ) is displayed for the relativistic 
FSUGold model of Ref. — a model that 
predicts an incompressibility coefficient for 
symmetric nuclear matter of K = 230 MeV. 
Note that the distribution of strength was 
obtained from a relativistic random-phase- 
approximation (RPA) approach as described 
in detail in Ref. (l2| . Further, the inset on 
Fig. CD shows a comparison of the theoretical 
predictions against the experimental centroid 
energies reported in Ref. [Hj]. Finally, the 
solid line in the inset provides a fit to the 
mass dependence of the theoretical predic- 
tions that yields E GMR (A) « [69 /A - 3 } MeV. 

The isoscalar monopole strength displayed 
in Fig.[T]is extracted from the low momentum 
transfer behavior of the longitudinal response 
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FIG. 1: (color online) Distribution of isoscalar 
monopole strength predicted by the FSUGold 
model of Ref. 10]. The inset includes a compar- 



ison against the experimental centroid energies 
reported in Ref. [l3], with the solid line provid- 
ing the best fit to the theoretical predictions. 

defined as follows: 



ft 1 (q,w) = X;k*»|p(q3l*o> S{ 



loj — u„ 



(1) 

Here \l/o is the exact nuclear ground state, 
\l/ n is an excited state with excitation energy 
u n , and p(q) is the Fourier transform of the 
isoscalar baryon density. That is, 



p(q) 



d 3 re -^(r) 7 roV'(r) , (2) 



where ip(r) is an isodoublet nucleon field, 7 
is the timelike (or zeroth) component of the 
Dirac gamma matrices, and To = 1 is the iden- 
tity matrix in isospin space. 

The important realization that the distri- 
bution of monopole strength in heavy nu- 
clei is sensitive to the density dependence 
of the symmetry energy has motivated a re- 
cent experimental study of the GMR along 
the isotopic chain in Tin. Indeed, the dis- 
tribution of isoscalar monopole strength in 
the neutron-even 112_124 Sn-isotopes has been 
measured at the Research Center for Nuclear 



Physics (RCNP) in Osaka, Japan [111, [14 
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This important experiment probes the in- 
compressibility of asymmetric nuclear mat- 
ter by measuring the distribution of isoscalar 
strength in a chain of isotopes with a neutron 
excess ranging from 6 = 0.11 (in 112 Sn) to 6 = 
0.19 (in 124 Sn). The experiment represents 
a hadronic compliment to the purely elec- 
troweak Parity Radius Experiment (PREX) 
at the Jefferson Laboratory that aims to mea- 
sure the neutron radius of 208 Pb accurately 
and model independently via parity- violating 
electron scattering (TH|, 16]. Such an accu- 



rate determination will have far-reaching im- 
plications in areas as diverse as nuclear struc- 



ture Il7j, heavy-ion collisions [la . Il9 
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231 ] . atomic parity violation 
and nuclear astrophysics j26|, 27, [2 

In Fig. [2] the experimental distribution of 
isoscalar monopole strength measured at the 
RCNP [ll|, 14] is compared against the pre- 
dictions of the hig hly successful NL3 (lol, HI 



and FSUGold [lOj models. As one is only in- 
terested in comparing the shape of the distri- 
bution and a particular ratio of its moments, 
the maximum of the theoretical curves - 
computed from the longitudinal response as 
described in the text — has been normalized 
to the experimental data. The A-dependence 
of the corresponding centroid energies is also 
displayed in Fig. [3] and compiled in Table [B 
Note that the centroid energy is computed 
from the ratio of the mi moment to that of 
the mo moment. That is, 



771 ""1 

"GMR = 



mi 
m 



L? S h (q ,uj)duj 



(3) 



where, consistent with the experimental anal- 
ysis (ill . [13 ]. the limits of integration have 
been chosen to be uji = 10 MeV and co> 2 = 
20 MeV. Further, to mimic the forward-angle 
experiment, the longitudinal response was 
evaluated at the "small" momentum transfer 
of g = 0.23 fin" 1 . 

A subtle telltale problem with Tin barely 
discernible in the inset on Fig. (TJ, be- 
comes magnified in Fig. [2] as one com- 
pares the experimentally extracted distri- 
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FIG. 2: (color online) Comparison between the 
distribution of isoscalar monopole strength in 
all neutron-even 112 Sn- Sn isotopes extracted 
from experiment (black solid squares) and the 
theoretical predictions of the FSUGold (blue 
solid line) and NL3 (green dashed line) models. 



bution of monopole strength against the 
theoretical predictions. While mean-field 
plus RPA calculations are typically un- 
able to describe the experimental width 
- which is in general composed of both 
an escape (particle-hole) and a spreading 
(multiparticle-multihole) width — such is not 
the case for the description of the centroid 
energies. Indeed, accurately calibrated mod- 
els, both non-relativistic [8j and relativistic 
(see Fig. CP, provide an adequate description 
of the GMR centroid energies in both 90 Zr 
(with 6 = 0.11) and 208 Pb (with 6 = 0.21) - 
nuclei with a neutron excess similar to those 
at the two extremes of the isotopic chain con- 
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continues to grow as the neutron excess in- 
creases? A stiff symmetry energy leads to a 
rapid softening of the nuclear incompressibil- 
ity This is the main reason behind the 
slightly larger (negative) slope displayed by 
NL3 relative to FSUGold in Fig. [3j The even 
larger (by more than 50%) slope displayed 
by the experimental data is unlikely to be 
solely related to the stiffness of the symme- 
try energy, as NL3 already predicts a neutron 
skin thickness in 208 Pb that appears overly 
large 
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FIG. 3: (color online) Comparison between the 
GMR centroid energies (mi /mo) in all neutron- 
even 112 Sn- Sn isotopes extracted from exper- 
iment (black solid squares) and the theoretical 
predictions of the FSUGold (blue up-triangles) 
and NL3 (green down-triangles) models. Also 
shown (filled red circle) is the result of Ref. [3] 
for the case of 116 Sn. 



Nucleus 


NL3 


FSUGold 


Experiment 


112 Sn 


16.98 


16.45 


16.2±0.1 


114 Sn 


16.92 


16.38 


16.1±0.1 


116 Sn 


16.81 


16.27 


15.8±0.1 


118 Sn 


16.70 


16.15 


15.8±0.1 


120 Sn 


16.66 


16.14 


15.7±0.1 


122 Sn 


16.54 


16.07 


15.4±0.1 


124 Sn 


16.43 


15.97 


15.3±0.1 



TABLE I: Giant Monopole Resonance centroid 
energies (in MeV) computed from the ratio of 
moments (mi /mo) as described in the text. All 
moments were obtained from integrating the dis- 
tribution of strength over the 10 < u < 20 MeV 
interval. 

sidered here. Why is then that both non- 
relativistic [ll|, 14, and relativistic mod- 
els consistently overestimate the centroid en- 
ergies in the Sn-isotopes? Or more colloqui- 
ally, why is Tin so soft ? And why is that the 
discrepancy between theory and experiment 



So why is Tin so soft and why does it 
become even softer with an increase in the 
neutron excess? Could there be a system- 
atic error in the experimental extraction? 
While possible, this is unlikely as an earlier 
independent measurement on 116 Sn [l3[ ap- 
pears to confirm the present (RCNP) result 
(see Fig. [3]). Could the GMR in Tin probe 
physics that has not been already constrained 
by nuclear observables? This also appears 
unlikely as existing density functionals are 
successful at describing a host of ground- 
state observables as well as collective exci- 
tations — including the GMR in 90 Zr, 144 Sm, 
and 208 Pb (see Fig. □ and Ref. 0). Could 
Tin be sensitive to pairing correlations and 
more complicated multiparticle-multihole ex- 
citations? The answer appears to be neg- 
ative [32], but even if it would be positive, 
why should Tin be sensitive to these effects 
but not Zr, Sm, and Pb? Clearly, the distri- 
bution of isoscalar monopole strength in the 
Sn-isotopes poses a serious theoretical chal- 
lenge, perhaps suitable for the new Univer- 
sal Nuclear Energy Density Functional (UN- 
EDF) initiative. Whatever the theoretical 
approach, however, one must remember that 
the challenge is not solely to describe the dis- 
tribution of monopole strength along the iso- 
topic chain in Tin, but rather to do so with- 
out sacrificing the enormous success already 
achieved in reproducing a host of ground- 
state properties and collective modes. 
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